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Abstract 

A Bethe-Salpeter-Faddeev (BSF) calculation is performed for the pentaquark + in 
the diquark picture of Jaffe and Wilczek in which Q + is a diquark-diquark-s three-body 
system. Nambu-Jona-Lasinio (NJL) model is used to calculate the lowest order diagrams 
in the two-body scatterings of sD and DD. With the use of coupling constants determined 
from the meson sector, we find that sD interaction is attractive in s-wave while DD 
interaction is repulsive in p-wave. With only the lowest three-body channel considered, 
we do not find a bound ^ + pentaquark state. Instead, a bound pentaquark + with ^ 
is obtained with a unphysically strong vector mesonic coupling constants. 
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1 Introduction 



The report of the observation of a very narrow peak in the K + n invariant mass distribution 
[1, 2] around 1540 MeV in 2003, a pentaquark predicted in a chiral soliton model [3], 
triggered considerable excitement in the hadronic physics community It has been labeled 
as + and included by the PDG in 2004 [4] under exotic baryons and rated with three 
stars. Very intensive research efforts, both theoretically and experimentally, ensued. 

On the experimental side, practically all studies conducted after the first sightings were 
confirmed by several other groups produced null results, casting doubt on the existence 
of the five-quark state [5, 6]. Subsequently, PDG in 2006 reduced the rating from three 
to one stars [4]. More recently, the ZEUS experiment at HERA [7] observed a signal for 
G+ in a high energy reaction, while HI [7], SPHINX [8] and CLAS [9] did not see it. This 
disagreement between the LEPS [1] and other experiments could possibly originate from 
their differences of experimental setups and kinematical conditions. So the experimental 
situation is presently not completely settled [10, 11, 12]. 

Many theoretical approaches have been employed, in addition to the chiral soliton 
model [3], including quark models [13], QCD sum rules [15], and lattice QCD [16] to un- 
derstand the properties and structure of G + . Several interesting ideas were also proposed 
on the pentaquark production mechanism. Review of the theoretical activities in the last 
couple of years can be found in Refs. [17, 18]. 

One of the most intriguing theoretical ideas suggested for + is the diquark picture of 
Jaffe and Wilczek (JW) [19] in which + is considered as a three-body system consisted 
of two scalar, isoscalar, color 3 diquarks (D's) and a strange antiquark (s). It is based, 
in part, on group theoretical consideration. It would hence be desirable to examine such 
a scheme from a more dynamical perspective. 

The idea of diquark is not new. It is a strongly correlated quark pair and has been 
advocated by a number of QCD theory groups since 60's [20, 21, 22]. It is known that 
diquark arises naturally from an effective quark theory in the low energy region, the 
Nambu-Jona-Lasinio (NJL) model [23, 24]. NJL model conveniently incorporates one 
of the most important features of QCD, namely, chiral symmetry and its spontaneously 
breaking which dictates the hadronic physics at low energy. Models based on NJL type of 
Lagrangians have been very successful in describing the low energy meson physics [25, 26]. 
Based on relativistic Faddeev equation the NJL model has also been applied to the baryon 
systems [27, 28]. It has been shown that, using the quark-diquark approximation, one 
can explain the nucleon static properties reasonably well [29, 30]. If one further take the 
static quark exchange kernel approximation, the Faddeev equation can be solved analyti- 
cally. The resulting forward parton distribution functions [31] successfully reproduce the 
qualitative features of the empirical valence quark distribution. The model has also been 
used to study the generalized parton distributions of the nucleon [32]. Consequently, we 
will employ NJL model to describe the dynamics of a diquark-diquark-antiquark system. 
To describe such a three-particle system, it is necessary to resort to Faddeev formalism. 

Since the NJL model is a covariant-field theoretical model, it is important to use 
relativistic equations to describe both the three-particle and its two-particle subsystems. 
To this end, we will adopt Bethe-Salpeter-Faddeev (BSF) equation [33] in our study. For 
practical purposes, Blankenbecler-Sugar (BbS) [34] reduction scheme will be followed to 
reduce the four-dimensional integral equation into three-dimensional ones. 

In Sec II, NJL model in flavor SU(3) will be introduced with focus on the diquark. The 
NJL model is then used to investigate the antiquark-diquark and diquark-diquark interac- 
tion with Bethe-Salpeter equation in Sec. III. In Sec. IV, we introduce the Bethe-Salpeter- 
Faddeev equation and solve it for the system of strange antiquark-diquark-diquark with 
the interaction obtained in Sec. III. Results and discussions are presented in Sec. V, and 
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we summarize in Sec. VI. 



2 SU(3) f NJL model and the diquark 

The flavor SU(3)f NJL Lagrangian takes the form 

C = $(i$-m)il) + Ci, (1) 

where ip T = (u,d,s) is the SU(3) quark field, and m = diag{m u ,m l j i ,m s ) is the current 
quark mass matrix. Cj is a chirally symmetric four-fermi contact interaction. By a Fierz 
transformation, we can rewrite Cj into a Fierz symmetric form £>i, q q = + ^(Cj)), 

where T stands for the Fierz rearrangement. It has the advantage that the direct and 
exchange terms give identical contribution. 

In the qq channel, the chiral invariant Ci >q q, is given by [35] 



<N 



Gi 



G 3 (^A» 2 + (^7 5 A /^ 



Go 



-G 4 



x q, 



+ 



(2) 



where a = ~ 8, and A° = yf/. If we define G 5 by -G 5 (^7%-) 2 = -(G 2 + G 3 + 
G 4 )('0i7 M AjV'i) 2 - G 2 (^7 M AjVi) 2 where i,j = u, d, then G 3 ,G 4 ,G 5 are related by G 5 = 
G 2 + with Gt, = G 3 + G 4 . In passing, we mention that the conventionally used G u 
and G n are related to G 5 , G v by G w = 2G 5 and G p = 2G 5 _- ±G V . 

For the diquark channel we rewrite into an form (%[}A'4) T ){'ip T Btp), where A and £> 
are totally antisymmetric matrices in Dirac, isospin and color indices. We will restrict 
ourselves to scalar, isoscalar diquark with color and flavor in 3 as considered in the JW 
model. The interaction Lagrangian for the scalar-isoscalar diquark channel [36, 37] is 
given by 

(3) 



d, s = G s ^C)X)^ T V T (C-V)A^ C A V> 



• n 2 • 

Z7 U 7 Z is 



where /3 A = y^X A (A = 2,5,7) corresponds to one of the color 3 C states. C 
the charge conjugation operator, and X's are the Gell-Mann matrices. 

The Bethe-Salpeter (BS) equation for the scalar diquark channel [36, 37] is given by 



—- I tr[(C- 1 7 5 T / 2 /3 A ))5(fc + q)(7 5 CT 2 f p A )S T (-q)}r s (q), 



(4) 



where the factors 4 and 2 arise from Wick contractions. S(k) = Qt — M + ie)^ 1 with 
M = M u = Md, the constituent quark mass of u and d quarks, generated by solving the 
gap equation. T s (q) is the reduced t-matrix which is related to the t-matrix by t s (q) = 
(7 5 Cr)/3^)r s (g)(C- 1 7 5 ^//? c A ). The solution to Eq. (4) is 



4iG, 



l + 2G s n s (<z 2 )' 



with 



IL((? 2 ) = 6i 



(2tt)' 



tr D tfS(qh 5 S(k + q)]. 



The gap equation for u, d and s quarks are given by 

Mi = mi- 8G1 < q~iqi >, 



(5) 



(6) 



(7) 
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with 



< <Mi >= -iN ( 



(8) 



where i = u,d, s. 

The loop integrals in Eqs. (6) and (8) diverge and we need to regularize the four- 
momentum integral by adopting some cutoff scheme. With regularization, we can solve 
the gap equation and t-matrix of the diquark in Eqs. (5) and (8) to determine the 
constituent quark and diquark masses. However, since our purpose in this work is not an 
exact quantitative analysis but rather a qualitatively study of the interactions inside G + , 
we will not adopt any regularization scheme and simply use the empirical values of the 
constituent quark masses M = M u ^ = 400 MeV, M s = 600 MeV, and the diquark mass 
M D = 600 MeV as obtained in the study of the nucleon properties [27, 28, 29, 31, 32]. 



3 Two-body interactions for strange antiquark- 
diquark (sD) and diquark-diquark (DD) channels 



In the JW model for G + , the two scalar-isoscalar, color 3 diquarks must be in a color 
3 in order to combine with s into a color singlet. Since 3 is the antisymmetric part of 
3 x 3 = 3 © 6, the diquark-diquark wave function must be antisymmetric with respect to 
the rest of its labels. For two identical scalar-isoscalar diquarks [ud]o, only spatial labels 
remain so that the spatial wave function must be antisymmetric under space exchange and 
the lowest possible state is p-state. Since in JW's scheme, G + has the quantum number 
of J p = i + , s would be in relative s-wave to the DD pair. Accordingly, we will consider 
only the configurations where sD and DD are in relative s- and p- waves, respectively. 

We will employ Bethe-Salpeter-Faddeev equation [33] to describe such a three-particle 
system of sDD. For consistency, we will use Bethe-Salpeter equation to describe two- 
particles subsystems like sD and DD, which reads as, 



where B is the sum of all two-body irreducible diagrams and Go is the free two-body 
propagator. In momentum space, the resulting Bethe-Salpeter equation can be written as 



where Go is the free two-particle propagator in the intermediate states, k and P are, 
respectively, the relative and total momentum of the system. 

In practical applications, B is commonly approximated by the lowest order diagrams 
prescribed by the model Lagrangian and will be denoted by V hereafter. In addition, it 
is often to further reduce the dimensionality of the integral equation (10) from four to 
three, while preserving the relativistic two-particle unitarity cut in the physical region. 
It is well known (for example, Ref. [38]) that such a procedure is rather arbitrary and 
we will adopt, in this work, the widely employed Blankenbecler-Sugar (BbS) reduction 
scheme [34] which, for the case of two spinless particles, amounts to replacing Go in Eq. 



T = B + BG T, 



(9) 



T(k',k;P)=B(k',k;P)+ / d 4 k" B(k' ,k" '; P)G (k" ; P)T(k" ,k; P), 



(10) 



(10) by 



G (k,P) 



1 



1 



(P/2 + k) 2 



— m 



\ {P/2 - k) 2 - ml 
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x «5<+> ((P'/2 + A;) 2 - m 2 ) <5« ((P'/2 - /c) 2 - mi) 



with 



G B bS( r H x = Si(|fc|) + fl2(|fc|) : 

2£ 1 (|£|)i<: 2 (|fc|) a-C^dfeD+^ad^ + ie' 



where s = P 2 and P' = yV/sP. The superscript (+) associated with the delta func- 
tions mean that only the positive energy part is kept in the propagator, and J5i2(|fc|) = 

k 2 + m 2 2 • 



3.1 sD potential and the t-matrix 

In Fig. 1 we show the lowest order diagram, i.e., first order in £i t qq in sD scattering. Due 
to the trace properties for Dirac matrices, only the scalar-isovector (ipX'jip) 2 , the vector- 
isoscalar (^7^ X^ip) 2 , and the vector-isovector (^7^ XJip) 2 will contribute to the vertex 
r. Furthermore, the isovector vertex (ipT\ a rip) 2 will not contribute since the trace in 



P-si r=Xf or Af Y (a=o,s 

-7 



< • i s 




Figure 1: sD potential of the lowest order in Ci t qq. 

^8 f\a\j.„ J\2 \o\2- 



flavor space vanishes, X)a=o(^/)33^ r /(^/^/^/) = Thus only the vector-isoscalar term, 
('tp^f^X^tp) 2 , remains. 

For the on-shell diquarks, the lower part of Fig. 1 which corresponds to the scalar 
diquark form factor, can be calculated as 

(p^+PDfTF^q 2 ) = i J ^ tr [( 5D C'-VA^ c A )5(A ; + g)7^(A ; )( 5D7 5 C'A 2 /3 c A )5 T (A;- m )] 

2 f d A k 



= Sigh J T^vtM^ + qh»S(k)S(pix - fc)], (13) 

where we have made use of the relations C _1 ( 7 M ) T C = —7^, tr c [f3^l3^} = 3. go is defined 
by 

with 



Hd(p 2 ) = 6i / ^j«r[5(fc)5(p - k)], (15) 
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and Mjj is the diquark mass. F v (0) is normalized as 2p^F v (0) = , such that 

F„(0) = 1. 1 

Then the matrix element of the potential VsD can be expressed as 

< SfDf\V\siDi > = (-v(psi))(-iV S D)(PDi,PDf)v(Psf) 

= (+m)(-G v )(-v(p si )h tl v(psf) [(Aj)33 • tTf (tf(tf) 2 )' 

x (P^+P^r^^-y (16) 

i.e., 

VsD = ^G v F v (q 2 )V- sD (p DhPDf ), (17) 

with 

VsDiPDuPDf) = ty Di +P' D f)/2- (18) 

Here the factor +16i in Eq. (16) arises from the Wick contractions, and the factor 
trf((\ 2 j) 2 ) in Eq. (16) is introduced to divide F v {q 2 ), since the factor trf((X 2 ) 2 ) is already 
included in the expression of F v (q 2 ) by a trace in flavor SU(3)f space. 




Figure 2: The BS equation for sD. 
The three-dimensional scattering equation for the sD system is now given by 

tsD(PDi,PDf) = VsD(PDi,PDf) 

+ 47T y dW ^f \ d Xi G^ S (\p D \,S 2 )K gD (\PDi\, \pL\^i)tsD(PD,PDf), 

(19) 



where X; = p Di ■ p D , p = p/\p\, s 2 = {p D i + Psi) = (PDf + Psf) , P D i = VPDi + M D> 



P°Df = ^P% f + M 2 D <m& 

KgD(\PDi\,\p D \, x i) = ^GvF v ((p' D - PDi) 2 )Kg D (p Di ,p' D )\ _ < , 2 2 , 

O ' J D — V ' D ~ Tm D 

KsD(PDi,P D ) = (#Di+P'D){-p'-s +M a )/2, 

with M s being the constituent quark mass of s and s. 



1 In the actual calculation we use the dipole form factor, F v (q 2 ) = (1 — q 2 /A 2 ) 2 with A = 0.84 GeV since 
the q 2 dependence for F v (q 2 ) in the NJL model is not well reproduced. 
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We also present the results for the interactions between diquark and u or d, which 
would be of interest when we study non-strange pentaquarks. One can just repeat the 
derivations we describe in the above and easily obtain 

Vud = V dD = - 16G 1 F s (q 2 ) + 32G 5 F v (q 2 )Vs D (p Di ,p Df ), (20) 

in analogous to Eqs. (17) and (18). 

We add in passing that, within tree approximation, the sign of the potential for sD is 
opposite to that of V s d due to charge conjugation, i.e., 

V s D{PDf,PDi) = -V sD {pDi,VDf)- (21) 
We can immediately write down the scattering equation for the sD as, 

tsD(pDf,PDi) = V sD (p Df ,p Di ) 

+ 4vr| dW ^f \f^ dx f G^ s (\p D \,s 2 )KM\PDf\Ap D \,x f )tMPD,PDi), 

(22) 

where x f =p Df -p D , Gf% s (\p D \, s 2 ) = G?% s (\p D \, s 2 ), and 

KsD(\pDf\,\PD\,Xf) = ^G v F v ((p' D - PDf) 2 )K sD (p D f,p' D )\ p ,^ = ^ ^ 2+M2 ^ 

K S D(PDf,p' D ) = -^ Df +p- D Ws+M s )/2, (23) 
with p' s = p'- s . 



3.2 Representation in /o-spin notation 

In the sD (or sD) center of mass system the wave function which describes the relative 
motion in J = ^, is given by the Dirac spinor of the following form (see [39, 40]), 



- (4ifS,i)K (24) 

\D „ (n° t?\ - ( V-Ps<t>s2{p%\Ps\) \ 

*8d(p°„Ps) = ^ sD (p° s ,Psh°, (26) 
VsD(P°s,Ps) = ^ sD (plPsh°, (27) 

where po = -p a = -ps, i.e., $> s d(p®,Ps) = ^ s d(p° s , ~Pd) and ^ sD (p%p s ) = ^sd(p®, ~Pd)- 
In the following we simply write p'q = IPqI^q^/) = \PQi(f)\ , Q = s,s or D. Note that 
the index 1 (2) corresponds to large (small) components for both s and s quark spinors. 

For a discretization in spinor space, we define the complete set of p-spin notation 
([39, 41]) for the operators O s d = V S D,t S D, V s d and K, s d = K s d, K s o of sD: 

O sD ,nm{PDf,PDi) = tr[ni(p s f)O sD (p D f,p Di )n m (p s , i )], (28) 
1C S D,nm(PDf,p'Di x f) = tr [^i(Psf)^sD{PDf,PD: x f)^m{p s )], (29) 
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where n,m= 1,2, fii(p) = and ^(p) = 7"P^> ^ = ^i(p) an d ^(p) satisfy 

tr[nU P )n m (p')} 

— Onl^ml + P ' P 3n2$m2- 

Concerning the sD spinor, the large and small components can be reversed by 7 s , 
with the minus sign which comes from the definitions Eqs. (25) and (27): ^> s dO^> s d = 
—^ S Dj 5 Oj 5 ^sD- Then we can define /9-spin notation for sD i.e., O s d = VsD,tsD,VgD 
and JCsD = KsD,K§d, 

OsD,nm(PDi,PDf) = -tr[^ n {psi)T '0 S D(PDi,P D fh^miPsf)], (30) 
IC§D,nm(PDi,PDi x i) = ~ tr [^niPsijl^K-sD (PDi, Pd, Xi)^VL m (p'^)] . (31) 

From Eqs. (19,22,28-31), each component n (n = 1,2) of spinors for the sD satisfy 
the following quadratic equation: 



0sn{Psi)tsD,nm(PDuPDf)(l>sm{Psf) = 4>tn(Psi) V§D,nm{PDi, PDf) 
(2vr 



2 f dp' i 1 f l i 

+47T^ J J2^PE>2 J i dx i G 'sD S '(PDi S2)KsD,nl(PDi,PD, x i) t sDM(PD,PDf)\(psm(Psf)- 



(32) 



A similar equation can be obtained for the sD by exchanging i <-> / and s <-> s in Eq. 
( 32 )- 

The explicit expressions of the p-spin notation for V^s)d and K g ^ D are given in 
appendix B. We note that there are important relations: 

VsD,nm(P,q) = -V s D,nm(P,q), 

VsD(p,q) = -V sD (p,q), 

KsD,nm(\p\,\<l\i x pq) — KsD,nm(\p | , \q | , Xpq), 

KsD{\p\,\q\,x pq ) = -K sD {\p\,\q\,Xpq). 

By the partial wave expansion in Eq. (69) in appendix A, the BS equation for t S D,nm 
in Eq. (32) for s-wave can be written as 

t l IS^n(PDi,PDf) = V*h D =°(PDi,PDf)+^ J T^dE^^^ 

(33) 



3.3 DD potential and t-matrix 

In the case of DD interaction, the lowest order diagrams are depicted in Figs. 3(a) and (b), 
with (a) the quark rearrangement diagram and (b) of the first order in £i iq q, respectively. 

We first show that the quark exchange diagram in Fig. 3(a) does not contribute due 
to its color structure, where a ~ d and i ~ I denote the color indices of the diaquarks and 
quarks, respectively. Since each diquark is in the color 3 [19, 36], the color factor for the 
qqD vertex is proportional to e a ij . Hence the color factor of the quark exchange diagram 
is given by 

taijtbiktclktdlj = Sabred + 5ad&bc- (34) 

As we discussed earlier, the color of the DD pair inside + is of 3 in order to combine 
with s to form a color singlet pentaquark. As color 3 state is antisymmetric under the 
exchange between diquarks in the initial and final states, the matrix element of Eq. (34) 
vanishes. 
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Figure 3: Lowest order diagrams in DD scattering. 



For the contact interaction diagram Fig. 3(b), only the direct term is shown since 
the exchange term does not contribute as it has the same color structure as the quark 
rearrangement diagram of Fig 3(a). It is easy to see that the color structure of Fig. 3(b) 
is proportional to S a b5 cc i. Then the terms in the interaction Lagrangian in Eq. (2) that 
can give rise to non-vanishing contributions are: 

GiOM^) 2 , -GatyfA^) 2 , -G V {^\ Q ^)\ (35) 

with a = ~ 8. 

We next calculate the form factors, which diagrammatically correspond to the lower 
part of diagram in Fig. 1. For V = 7^ Ay, we obtain 



tr f (\ a f {\)f){pm+VD f ) 

I I I 1 I 



„ F v {q 2 ) 



2, /1„ , ,„„ , 2 



3^0 + ^3^81 (PDi+PD f rF v (q z ), (36) 



and for T = A^, we get 

*-/(w) 5 ^)-(yi*.+ 1 /i«-)«tf 



(37) 



where the factor tr/((A^) ) in Eqs. (36) and (37) is introduced by the same reason for 
Eq. (16), and we have used tr(X 2 f X a f X 2 f ) = 2(y / f<5 a0 + ^ 



For the on-shell diquarks, F s {q ) is calculated as 



d A k 



W) = ij j^tr[(g D C- L 1 b Xjp A )S(k + q)S(k)(g D ^C\jp A )S 1 (k-p Di )] 



, / d A k 



= Ugh J j—utr[S(k + q)S(k)S(k- PDi )}. (38) 
With the form factors F v (q 2 ) and F s (q 2 ) obtained in the above, Vdd is given by 



-iV DD (p Di ,p D f) = +128i 
= 128* 



GlF 2 (q 2 ) - (g 2 + (pDli+PDlf) ■ (PD2i+PD2f)F 2 (q 2 ) 

,-i2/„2\ /~i /„_ i _ \ /„_ i „ \T?2r2 



G 1 F s \q 2 ) - G 5 ( PDli +p Dlf ) ■ ( P D2 l +PD2f)F v 2 (q 2 ) , (39) 



2 Same as the case for sD potential, we use the dipole form factor, F s (q 2 ) = c s (l — q 2 /A 2 ) 2 with A = 0.84 
GeV and c s is a constant. In the original NJL model calculation with the Pauli-Villars (PV) cutoff, c s is given 
by F s (0) = c s = 0.53 GeV [32]. 
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where the factor +128i in a first line of Eq. (39) comes from the Wick contractions, and 
in a second line we have used the relation between couplling constants in meson sectors; 
G5 = G 2 + \G V which is explained in section 2. The momenta of the diquarks in the 
initial and final states in Fig. 4 are given by 

VDli(f) = (v^A PDi(f)), 

Pmnf) = (>/*2 A -&*(/))> (40) 
with q = Pmf-PDli = PD2i~PD2f- s 2 = 4(p^ + M|,) = 4{p^ f + M^) is the DD center 
of mass energy squared. 



P Dli P Dif P/2+p' 




Figure 4: BS equation for DD. 

As in the case of sD scattering, we use the BbS three-dimensional reduction scheme 
and the resulting equation for DD scattering reads as 

f d ?J p' 1 i i 

tDD(PDf,PDi) = V DD (p Df ,p Di )+ I j^—^V DD (p Df ,p')G^ b S(\p'\, s 2)tDD{p' ,PDi), (41) 



with 



G B D b S(\p'\,S 2 ) : ' 



AE D {\p'\)( S2 /A-E D {\p'\Y + ie) 

±E D {\p'\)(p* f -p'* + iiy (42) 

with E D (\p'\) = yfiri + M*. 

In the JW model for + , the diquark-diaquark spatial wave function must be anti- 
symmetric and we will consider here only the lowest configuration, namely, DD are in 
relative p-wave. Partial wave expansion of Eq. (69) then gives 

dp' 
(2^f 



%h(pf,Pi) = Vh = D 1 (Pf,P l ) + ^p' 2 G^(p', s 2 )V}T D 1 (pf, P ')t l ^(p',p i ), (43) 
with p i(/) = \p Di ( f )\,p' = \p'\- 



4 Relativistic Faddeev equation 
4.1 3-body Lippmann-Schwinger equation 

For a system of three particles with momenta k^s (i = 1,2,3), we introduce the Jacobi 
momenta with particle 3 special choice: 

k\ = fi 1 P + p + ai q 3 
k 2 = fi 2 P -p + a 2 q 3 

h = ^P + a 3 q 3 , (44) 
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with Mn = 1 an d a 3 = —d\ — a 2 . For the coefficients we find [i n = m n /M, M = 

7711 +771,2 +7713, and CKl = 777i/t77i2 , a 2 = m 2 /mi 2 , Q3 = — 1, where TTljj = 777j + 777j {% ^ j). 

In terms of the Jacobi momenta the total kinetic energy is given by: 

p2 c: 2 ~2 

— + ^ + , (45) 

2M 277712 2777( 1 2)3 

where = rrikmij / M . 

New integration variables are chosen to be: p = f p3 p with f p3 = \J2m\ 2 and q 3 = f q3 q 
with f q3 = ^/2m (12 )3, and in general for cyclic (ij/c), / pi = ^277i jfc and f qi = ^2m {jk)i . 
In terms of the new integration variables we have 

Ktot = ^+p 2 +q 2 , (46) 
and the 3-body Lippmann-Schwinger equation for the T-matrix becomes: 

Tip, q) = V + / p3 3 J -0jsj -0js V G 3b', q') W, q% (47) 

with G 3 (p, q) = l/(z — K t ot)- The parameter z is implicit in the arguments of T and G 3 
in Eq. (47), a convention to be followed hereafter. 

Similarly we define the Jacobi momenta pi, qi with particle i as the special choice. The 
momenta are related to each other as 

Pi = dijpj + bijqj, qi = Cijpj + dijifj, (48) 



where (ijk) are cyclic, and = —[mirrij / {rrii +mf c )(mj + 777; s )] 1 / 2 , bij = ^Jl — afj = —bji, 
It can be shown that the total angular momentum is related to the angular momentum 

'pi 

and l q i by 

3 3 

L = 2 ( f i x = Yl (> + + £■ (49) 
i=i i=i 

With these three choices of Jacobi momenta we may introduce corresponding 3-particle 
states I > n where particle n plays a special role. For the 3-particle T-matrix we have 

< h, k 2 , k 3 \T\a >=„< p n , q n \T\a >, (50) 

or in terms of the Faddeev amplitudes T n , 

< k 1} k 2 , k 3 \T\a > = Ti(pi,gi) + T 2 (p 2 , q 2 ) + T 3 (p 3 , q 3 ), (51) 

with T n (p n ,q n ) = n < p n ,q n \T n \a >. 

For the pentaquark system we now chose particles 1 and 3 as the diquark and particle 
2 to be the s. The Faddeev equations for T = T x +T 2 + T 3 with T i = t i + Y J tiG 2 (s)Tj (i = 

1,2,3), with ti denoting the two-body t-matrix between particle pair (jk), become 
Ti(pi, qi) = f P 3 3 f q3 3 j ^13 G 3 (p' 3 ,q' 3 ) T 3 (p 3 ',q 3 ') 

+ fp^f^J §$J^$ K ^ G s(pU) T 2 (p 2 ',q 2 '), (52) 

where the channels 1 and 3 correspond to D(sD) states and channel 2 to the s(DD) 
states. Since diquarks obey Bose-Einstein statistics, we have T 3 (p 3 ,q 3 ) = Ti(—p 3 ,q 3 ) and 
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T 3 (p 3 ,q 3 ) = Ti(—pi,qi). We note that the symmetry property which requires the ampli- 
tude T be anti-symmetric with respect to interchange of the 2 diquarks is automatically 
satisfied by the angular momentum content L = l q ± = l p2 = 1, l p i = l q2 = 0. 
The s(DD) T-matrix T 2 satisfies 

d 3 p'i f d 3 Qi 
(2^)3 J J27 

The kernels Ki 3 and K\ 2 are expressed in terms of the sD t-matrix 

(2vr) 3 



T 2 (p 2 ,q 2 ) = 2fpffqf J ^ J |^ K 21 G^q',) Ti(pi',qi 



(53) 



K 13 = K 12 = tg D (pi,Pi '\z-q-y) 



f 3 
Jqi 



Similarly the kernel K 2 \ is given by 



K 2X = t DD (p 2 ,p 2 '; z-ql) ^-S^[q 2 



f 3 
Jqi 



1 - 91 



92 



The term with K\ 3 can be worked out by making use of the 5-function relation 

(qi ~ Qi 2 ) $ (cos 9 q3 - cos e q i^j 5 (<^ - 4> q3 



5 (3) [qi-qi'] 



9i 



(54) 



(55) 



(56) 



and the linear relation q\ ' = c\ 3 p 3 ' + ^139*3 ', which lead to 



5^[q,-qi'] = 



1 



91013^13^393 



-5 cos 9„i „i — 



12 _ 12 12 _ j2 12 
9l C 13P3 "1393 



p 3«3 2c 13 d 13 p' 3 q' 3 



5 (cos 9 q3 - cosflg/) 5 (<^> ? / - ^ 3 ) . 



We mention that similar expression for a delta function in the term Ki 2 can also be 
obtained by replacing 3^2. 

Performing a partial wave expansion for the D(sD) amplitude 

Ti(pi,q\) = 4nYf pi0 (n pi )Y lqi0 (n qi )T l L (p 1 ,q 1 ), (58) 

and for the sD t-matrix t S D(pi,Pi '] z — q\), 

tsoipup! '■ z-q\)= 4nYr piO (n pi )Y lpiO (n p 0t% 1 \p 1 ,p\-,z - qf), (59) 

yield 

Tf(pi,qi) 

= c 3 / q' 3 dq> 3 / Ps'dp's t% 1 \p 1 ,p' 1 ;z - qf) X 13 — — G 3 (p' 3 , q' 3 ) T 3 L (p> 3 , q> 3 ) 

Jo JA l3 c\ 3 di 3 qip 3 q 3 

' q 2 dq' 2 [ Bl2 p' 2 2 dp' 2 t^ 1 \ Pl ,p' 1 ;z- q 2 1 ) X 12 1 - - G 3 (p' 2 ,q' 2 ) T 2 L (p' 2 ,q' 2 ), 

JA 12 C12 «i2 91 p 2 q 2 

(60) 



-c 2 



with 



2 2 

C3 = —T=(fp3fq3/fql) , ^2 = —J={fp2fq2/fq\) , 
V 71 " V 71 " 



(61) 



and where the boundaries A, B for the p' integration can easily be found from the condition 
q\ = q[ 2 in Eq. (57), given by 



(62) 
(63) 



Aij — 


(HjQj + 9i 


dij 


Bij = 


Cijq'j ~ qi 
dij 
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For the s(DD) amplitude T 2 , partical wave expansion gives, 



T 2 (P2,Q2) = 2ci / q[ dq[ / p x dp\ 

JO JAoi 



x t^\p 2 ,p' 2 ;z- q\) X 21 1 - - G 3 (pi,9i) ?Wi,<Zi), (64) 

C21 "21 ?2 Pi <?i 

where A21 and -B21 are given by Eq. (63), and 

ci = -^(UiU/Uf. (65) 

In the above equations are angular momentum functions depending on the states 
we consider. In our case, the sD 2-body channel is a s-wave, Ip = 0, and the DD channel 
a p-wave, Ip = 1. Hence, for the 3-body channel with total angular momentum L = 1 
we have for the D(sD) 3-body channnel lp\ = 0, Iqi = L and Zp3 = 0, ^3 = L, while for 
s(DD) lp 2 = 1, ^2 = 0. The obtained JQj have the form 

x n = -^—^ Y iq3o(8q :i qi), X 12 = — —Y lq2 o(6 q2 qi ), X 21 = j—^=Yi P20 (e p2Pl ). (66) 
4.2 Relativist ic Faddeev equations 

Following Amazadeh and Tjon [42] (see also [33]) we adopt the relativistic quasi-potential 
prescription based on a dispersion relation in the 2-particle subsystem. Then the 3-body 
Bethe-Salpeter-Faddeev equations have essentially the same form as the non relativistic 
version. Taking the representation with particle 3 as special choice we may write down for 
the 3-particle Green function a dispersion relation of the (l,2)-system, i.e., 



E 1 (k 1 ) + E 2 (k 2 ) 
Ei(h) E 2 (k 2 ) s 3 -q 2 3 - {E 1 (k 1 ) + E 2 {k 2 )f 



G 3 {P3, q 3 ; s 3 ) = 171 TT„ \ 17 7Z ^ ~ J2 7 F 7TT\ I ip 7Z ^^2 ' ( 67 ) 



with Ei(k\) = ykf +mf, E 2 (k 2 ) = \jk\ + m 2 , and S3 = P 2 being the invariant 3-particle 
energy square. In the 3-particle cm-system we have ^/ss = M + E\>. The resulting 2-body 
Green function with invariant 2-body energy square s 2 has then the form of the BSLT 
quasi-potential Green function 

E^ki) + E 2 {k 2 ) 1 
° 2{P3] S2 > ~ E l (k l )E 2 (k 2 ) s 2 -(E l (k l ) + E 2 {k 2 )f {68} 

This quasi-potential prescription for G 3 has obviously the advantage that the 2-body 
t-matrix in the Faddeev kernel satisfies the same equation as the one in the 2-particle 
Hilbert space with only a shift in the invariant 2-body energy. So the structure of the 
resulting 3-body equations are the same as in the non relativistic case. 



5 Results and discussions 

In the NJL model some cutoff scheme must be adopted since the NJL model is non- 
renormalizable. However, in this work we will not use any cutoff scheme but simply 
employ the dipole form factors for the scalar and vector vertices. Namely, the NJL model 
is only used to study the Dirac, flavor and color structure of the sD and DD potentials. 

For the values of the masses M u ^, M s and Mo, we use the empirical values M = M u = 
M d = 400 MeV and M s = M D = 600 MeV [32]. We will treat the coupling constants 
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Gj (i = 1 ~ 5) in Eq. (2) as free parameters. For the sD channel, it depends only on 
G V = G 3 + G A = §(G 5 - G 2 ) as seen in Eq. (16). 

In the NJL model calculation with the Pauli-Villars (PV) cutoff regularization [32], 
the coupling constants G n , G p and G^ are related with the parameters used in our work 
by G\ = G-x/2, G2 = Gp/2 and G5 = G w /2. Thus by using the values of mesonic coupling 
constants in the NJL model, G v is determined as G v = §(G w /2 — G p /2) = | (7.34/2 — 
8.38/2) = —0.78 GeV -2 . We remark that the sign of G v is definitely negative since 
experimentally omega meson is heavier than the rho meson. Then the interaction between 
s and diquark in s-wave is attractive, as can be seen from the sD s-wave phaseshift shown 
in Fig. 5 with G v = —0.78 GeV -2 , while the interaction between s and diquark is repulsive 
which can be seen in Fig. 6. In both figures we find that the magnitudes of the phaseshift 
is within 10 degrees, that is, G v = —0.78 GeV -2 gives very weak interaction between s (s) 
and diquark. As we can see in Figs. 5 and 6, generally the phaseshift in s-wave is more 
sensitive to three momentum than that in p-wave. We note that sD and sD phaseshift are 
not symmetric around the pe axis, which can be understood from the decompositions of 
t S D and tsD i n the spinor space in appendix B. We further mention that if G v is determined 
from the A hyperon mass M\ = 1116 MeV within the sD picture, one obtains G v = 6.44 
GeV -2 , which is different from G v = —0.78 GeV -2 determined from meson sector in the 
NJL model in sign. In this case the rho meson mass is larger than the omega meson mass, 
that is, the vector meson masses are not correctly reproduced. 

DD phaseshift is plotted in Fig. 7 where we have used the values of coupling constants 
Gi = G^/2 = 5.21 GeV -2 and G 5 = G w /2 = 3.67 GeV -2 which are determined from 
meson sectors in the NJL model calculation with the Pauli-Villars cutoff [32]. We can 
easily see that the phaseshift <5; is definitely negative i.e., the DD interaction is repulsive, 
and its dependence on three momentum pe is very strong and almost proportional to pe 
both for s-wave and p-wave. This strong pe dependence of phaseshift comes from the p\ 
dependence of a second term (pou + PDif) ■ (PD2i +PD2f) in Eq. (39). 

The G v dependence of the sD binding energy, E s n, is presented in Fig. 8. We find that 
the sD bound state begins to appear around G v = —5 ~ —6 GeV -2 , becomes more deeply 
bound as G v becomes more negative. It is easily seen that E s d is almost proportional 
to G v . However even for the case of a weakly bound state with \E s d\ less than 0.1 GeV, 
it will require a value of — G v = 5 ~ 6 GeV -2 which is about eight times larger than 
the — G v determined from meson sector in the original NJL model with the PV cutoff 
regularization. 

For the calculation of the pentaquark binding energy we use the relativistic three-body 
Faddeev equation which is introduced in section 4. If the pentaquark state is in J p = \ 
state with which we are concerned in the present paper, the total force is attactive but 
there is no pentaquark bound state. 

On the other hand if the pentaquark state is in J p = \ state, a bound pentaquark 
state begins to appear when G v becomes more negative than —8.0 GeV -2 , a value in- 
consistent with what is required to predict a bound A hyperon with M\ = 1116 MeV in 
a quark-diquark model as mentioned in Sec. 5. The lowest configuration which would 
correspond to a J p = \ state is for the spectator s to be in p— wave w.r.t. to a DD 
pair in p— wave, or alternatively speaking, the spectator diquark in relative s-wave to sD 
in s-wave. Our results for the binding energy of a J p = \ pentaquark state for the 
case with and without DD channel are given in Table 1. It is found that although the 
DD interaction is repulsive, including the DD channel gives an additional binding energy 
which is leading to the more deeply pentaquark boundstate. It is because the coupling 
to the DD channel is attractive due to the sign of the effective kernel K21 in Eqs. (53, 
55). This depends on the recoupling coefficients X21, X\2 in Eq. (66) and the 2-body 
t-matrices. 
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Lr v [LreV \ 


h B {bq)[MeV \ 




-8.0 


47 


77 


-9.0 


87 


139 


-10.0 


132 


205 


-12.0 


226 


333 


-14.0 


316 


505 



Table 1: The binding energy of J p = | pentaquark state. E B (5q) (E B (5q)) is the binding 
energy without (including) the -D.D channel. 



In Fig. 9 (10) the phaseshift of sD is plotted, where the coupling constant is fixed 
at G v = -8.0 GeV~ 2 (G v = -14.0 GeV -2 ). It is easily seen that in Figs. 9 and 10 the 
phaseshift of sD in s-wave is positive for small pe < 0.3 GeV and pe < 0.45 GeV, but 
it changes the sign around pe = 0.3 and pe = 0.45 GeV, thus the phaseshift of sD in 
s-wave is very sensitive to three momentum pe- Whereas the phaseshift of sD in p-wave 
is definitely positive. 

In Fig. 11 we plot the phaseshift of sD with the coupling constant G v = —14.0 
GeV -2 which is same as the one used in Fig. 10. Different from the phaseshift of sD the 
phaseshifts of sD in s and p-wave do not change the sign for higher three momentum pe, 
i.e., the sign of the phaseshifts are definitely negative. 

From the above results we find that even if we use a very strong coupling constant G v 
which is unphysical because it gives much larger mass difference of rho and omega mesons 
than the experimental value, M u —M p = 13 MeV, it is impossible to obtain the pentaquark 
bound state with J p = i + . With only the J = \ three-body channels considered, we do 
not find a bound J p = ^ + pentaquark state. The J p = ^ channel is more attractive, 
resulting in a bound pentaquark state in this channel, but for unphysically large values 
of vector mesonic coupling constants. 

6 Summary 

In this work, we have presented a Bethe-Salpeter-Faddeev (BSF) calculation for the pen- 
taquark B + in the diquark picture of Jaffe and Wilczek in which B + is treated as a 
diquark-diquark-s three-body system. The Blankenbecler-Sugar reduction scheme is used 
to reduce the four-dimensional integral equation into three-dimensional ones. The two- 
body diquark-diquark and diquark-s interactions are obtained from the lowest order dia- 
grams prescribed by the Nambu-Jona-Lasinio (NJL) model. The coupling constants in the 
NJL model as determined from the meson sector are used. We find that sD interaction 
is attractive in s-wave while DD interaction is repulsive in p-wave. Within the truncated 
configuration where DD and sD are restricted to p- and s-waves, respectively, we do not 
find any bound \ + pentaquark state, even if we turn off the repulsive DD interaction. 
It indicates that the attractive sD interaction is not strong enough to support a bound 
DDs system with J p = ^ + . 

However, a bound pentaquark with J p = ^ begins to appear if we change the vector 
mesonic coupling constant G v from —0.78 GeV -2 , as determined from the mesonic sector, 
to around G v = — 8 GeV -2 . And it becomes more deeply bound as G v becomes more 
negative. 
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Figure 5: Three momentum pe dependence of the phaseshift 81 for the sD interaction with the 
coupling constant G v = —0.78 GeV~ 2 . 
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Figure 9: Three momentum pe dependence of the phaseshift 5i for the sD interaction with the 
coupling constant G v = —8.0 GeV -2 . 




Figure 10: Three momentum dependence of the phaseshift 5i for the sD interaction with 
the coupling constant G v = —14.0 GeV~ 2 . 
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Appendices 

A Partial wave expansion 

In the 2-body center of mass frame the partial wave expansion is defined by 

21 + 1 

tip f, Pi) = 4?r Pl(cos9 PlPf ) < p f l\t\pil > 

= ^(21 + l)Pi(cos9 PiPf )t l i\p f \, \m\), 
i 



(69) 



with Pi(f) =Pu(f) = ~P2i(f)- Then t i\pf\, \pi\) in Eq. (69) is written in terms of t(pf,pi) 

by 

i f 1 

t l i\Pfl\Pi\) = g J dcosO PiPf Piicos9 PiPf )t(pf,pi). (70) 



The phase shift 5i is given by 



t ip,p) = — e l sin5i, 

p 



(71) 



where p = \pu\ = \p 2 i\ = \pif \ = \p2f \ and s 2 = (pu +P2i) 2 = (pif +P2f) 2 - 

B The results for Vs( s )D,nm and K- s{s)D ^ nm (n, m = 1, 2) 

In this appendix we show the results for Vs(s)D,nm an d K g ^ D ^ nm (n, m = 1,2) defined in 
Eqs. (28-31): 



VsD,ll(PDi,PDf, 



X = 



Ppj+Ppf 
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PDf + XPDi _ Psf + XPsi 



2 2 

p Di + Xp Df _ p S i + XPsf 



VsD,12{pDi,PDf,x) = 
VsD,2l(pDi,PDf,x) = 

X 

2 



VgD,22(PDi,PDf,- 



-XPm+PDf), 



and 



K S D,ll(PDi,PDi x i) = 2 ^ + (\/«2 -Pd)(PJ0 + PDi) + PD + XiP D PDi 



K S D,12{PDi,p' D i x i) 
KsD,2l(PDi,PDi x i) 
KsD,22{PDhPDi x i) 



{p' D + XiPDi)(M s - + p£) - + Pd) 

(PDi + Xip' D )(M s + - p£) + Xip' D (p° Di + p£) 

XiM s (p° Di + p2) - (PDiP'd + «»Pjo) + x *(Pd - Vs2~)(P°Di + Pd) 



where x = p Di ■ p D f , Xi = p Di ■ p D . 
VsD,nm and K sD 

,nm are related with V s D,nm 

and K sDi nm by 



V S D ;nrn(Pi Qi Xpq) 
KsD,nmiPi 0.1 Xpq) 



V S D ,nm(.P ■> Qi Xpq\ 
K-sD,nm(Pi Qi Xpq). 



C Parametrizations for tsD an d t s o 

tsD can be parametrized as 

tsDipDhPDf) = A P [ F S P ' + F T P ' ia ^PDfPDi 

p,p'=± 

where A± = ^r 1 . Components of t S D is written as 



tsD(PDi,PDf) 



F s + + F T +i ° ■ " 



Fo~ + F^~ia ■ n 



(72) 



(73) 



where n = p£>/ x PDi,v = p D fPDi — PdiPdj, and ± means upper and lower components 
in the spinor space i.e., (tsD)p,p> = ^-ptsD^-p 1 - 

The decomposition into upper and lower components in eq. (30) for t S D gives 

tsD,ll(PDi,PDf) = ~ F s~, 

tsD,12(PDi,PDf) = -F^ + (xp Df PDi- P°DiPDf), 

tsD,2l(PDi,PDf) = -Fr~(PDfPDi - XPDiPDf), 

tsD,22(PDi,PDf) = -F^ + p D iPDf(x 2 ~ I). 

We can parametrize t s D in the same way (= eq. (72)) 



tsD(pDf,PDi) = A P F S P + F T ' i(J pvP D fPDi 

p,p'=± 



A 



■pi 



(74) 



where A± = ^p 1 . 



Similar to tgD the decomposition into upper and lower components by eq. (28) gives 
tsD,ll(PDf,PDi) = F^ + , 

tsD,12(PDf,PDi) = F^(p Df PDi-Xp°DiPDf), 
tsD,2l(PDf,PDi) = F^ + (xp Df pDi-p D iPDf), 
tsD,22(PDf,PDi) = F T ~p D iPDf(x 2 ~ I). 
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